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Abstract. 

We consider the possibility that the primordial curvature perturbation was generated 
through the curvaton mechanism from a scalar field with an electric charge, or precisely the 
Standard Model U(l) weak hypercharge. This links the dynamics of the very early universe 
concretely to the Standard Model of particle physics, and because the coupling strength 
is known, it reduces the number of free parameters in the curvaton model. The gauge 
coupling also introduces several new physical effects. Charge fluctuations are generated 
during inflation, but they are screened by electron-positron pairs therefore do not violate 
observational constraints. After inflation, the curvaton interacts with thermal radiation which 
destroys the curvaton condensate and prevents the generation of curvature perturbations, 
unless the inflaton dynamics satisfy strong constraints. The curvaton also experiences a 
period of parametric resonance with the U(l) gauge field. Using the standard perturbative 
approach, we find that the model can generate the observed density perturbation for Hubble 
rate -ff* > 10* GeV and curvaton mass m > 10~^H^, but with a level of non-Gaussianity 
(Inl ^ 130) that violates observational constraints. However, previous studies have shown 
that the parametric resonance changes the predicted perturbations significantly, and therefore 
fully non-linear numerical field theory simulations are required. 
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1 Introduction 

The inflationary paradigm is commonly accepted as the origin of the primordial fluctuations 
that seeded structure in the Universe. Models of inflation are usually based on quantum 
field theory, with fields unrelated to the Standard Model of particle physics. However, some 
interaction between the infiaton and the Standard Model fields is required in order to reheat 
the Universe, which means that a single theoretical framework for both inflation and particle 
physics is ultimately required. 

If the inflaton field is responsible for generating the primordial curvature perturbations, 
its coupling to other fields has to be extremely weak in order to avoid radiative corrections 
that would spoil the flatness of the potential. This makes it difficult to couple the infiaton 
to Standard Model gauge fields, which all interact relatively strongly. A less constrained 
alternative is the curvaton model [1-4], in which the perturbations are generated by a separate 
scalar field known as the curvaton. The curvaton is light and subdominant during inflation 
and gains an isocurvature perturbation. After infiation has ended, the energy density of 
the curvaton grows relative to the background radiation. When the curvaton finally decays, 
its isocurvature perturbation is converted to an adiabatic perturbation that can seed the 
structure in the Universe. 

The purpose of this paper is to investigate whether the curvaton field could be charged 
under a Standard Model gauge group. This would have the attractive feature that the 
properties of these interactions are known, in contrast with typical curvaton models, which 
have so much freedom that it is hard to make definite predictions. Of the three Standard 
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Model gauge groups, the most promising is the U(l) weak hypercharge, because the SU(2) 
and SU(3) groups are confining and would give more complicated physics. In contrast, a U(l) 
charge would essentially mean that the curvaton has an electric charge, and it is relatively 
easy to investigate the various constraints that arise from this. 
To be specific, we consider the Lagrangian 



where a is the curvaton field, g' is the Standard Model U(l) gauge coupling, is the 
Standard Model U(l) gauge field and Ff^i^ is the corresponding field strength tensor. The 
coupling strength is known [7], g' « 0.36, and thus the number of free parameters is reduced.^ 
For concreteness, we assume that the curvaton carries one unit of hypercharge, Y = 1. The 
hypercharge has to be an integer to allow the curvaton to decay into Standard Model particles, 
and a higher value would not change our conclusions significantly. For Y = 2, the curvaton 
could have a Yukawa coupling to right-handed electrons, but for other values the only other 
renormalisable term is a bilinear coupling a^a^^^ to the Higgs field. For simplicity, we 
assume that this coupling is negligible. The opposite case would essentially correspond to 
the model discussed in Refs. [5, 6]. 

To determine whether such a model is viable, we apply both theoretical and observa- 
tional constraints, for example requiring a stable vacuum and the correct amplitude of the 
curvature perturbation In this model, three mechanisms affect the curvaton's dynamics 
after inflation: non-perturbative parametric resonance into photons, interactions with the 
thermal bath, and perturbative decay. We find that a perturbative estimate of the non- 
Gaussianity (/nl) appears to rule out the model at 95% confidence level, but this ignores 
the parametric resonance which is known to have a significant effect on perturbations [5] and 
therefore it cannot be relied on. Instead, non-linear field theory simulations are ultimately 
needed to determine the perturbations. 

In Section 2, we give a brief overview of dynamics of the model. In Section 3, we discuss 
the charge fluctuations produced during inflation and show that, because of electric screen- 
ing, they are compatible with observations. Then in Section 4, we give the conditions that 
the model must satisfy to have the qualitative behaviour of a curvaton model. These are 
conditions such as requiring a light curvaton during inflation and no false vacuum in the cur- 
vaton potential. In Section 5, we then expand the discussion of the dynamics after inflation, 
analytically calculating timescales for each process. In Section 6 we discuss the generation 
of curvature perturbations and production of transient cosmic string-like structures. We 
conclude in Section 7. 

2 Overview of the model 

The electrically charged curvaton model follows the standard curvaton scenario in many ways, 
but with more complicated dynamics and tighter constraints. We first review the standard 
curvaton scenario. There are two fields: the inflaton and the curvaton. During inflation, the 
inflaton dominates the energy density, but its perturbations are assumed to be negligible. 
The curvaton fleld is light (compared to the Hubble rate^ H^), and therefore it develops 
a nearly scale-invariant spectrum of fluctuations, in the same way as the inflaton. After 

^The coupling will run to larger values — we quote the weak scale value. 
^We use subscript * to represent values at the end of inflation. 
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inflation, it is assumed that the inflaton reheats to produce a thermal bath, which becomes 
the dominant form of energy in the Universe. 

As the Universe expands, H decreases, and when H < m, the curvaton begins to 
oscillate in its potential. Assuming a harmonic potential,'^ its equation of motion is 

a + 3H& + m^a = 0. (2.1) 

In a radiation-dominated background, H{t) = l/{2t) and a{t) = yj2H^t. The evolution of a 
is then given by 

21/4 r(^) (7. 



0.86 \-r cos(mt-37r/8). (2.2) 



The energy density of the curvaton field. 



p.«mV(t)2«0.74^^^, (2.3) 



decreases as a ^, whereas that of the radiation. 



decreases proportional to a~^. Thus, the relative fraction of curvaton r{t) grows. 

Eventually, the curvaton decays to Standard Model particles which have the equation 
of state of radiation — at this point the perturbations in the curvaton become adiabatic. 
Assuming that this happens instantaneously at time tdeo the amplitude of the curvature 
perturbation is given by [3] 

^ _ r-(tdcc) ^P^ 

C~ — ^ , (2.6 

3 Per 

This needs to agree with the observed amplitude (" ~ 10~^. In the standard curvaton model, 
in which Eq. (2.3) remains valid until the decay time, this becomes 

C~ — 5 da* ^ — . (2.7) 

With the same assumption, the non-Gaussianity of the perturbations is given by 

In the electrically charged curvaton model, the dynamics between the end of inflation 
and C becoming adiabatic are non-trivial. The large coupling g' means that photon-curvaton 

^Other potentials have been discussed in Refs. [8, 9]. 
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interactions play an important role, through a number of processes. A thermal bath contain- 
ing photons will interact with the condensate. This thermal bath could either be produced 
at the end of inflation by reheating of the inflaton, or later by curvaton decay products. 
The fate of the model strongly depends on whether any interactions with the thermal bath 
lead to full chemical equilibrium. If they do, then becomes adiabatic when that happens. 
Otherwise, the relative fraction of curvaton can continue to grow until the final decay. 

As we will see in Section 5.2, the large photon-curvaton coupling also leads to resonant 
production of photons as the condensate oscillates. This changes /Oo-, and makes it a more 
complicated function of a^:. Therefore Eqs. (2.7) and (2.8) are no longer valid. If the resonance 
ends before backreaction makes the dynamics non-linear, some condensate will remain, and in 
principle the dependence can be calculated from linear theory. Previous numerical simulations 
indicate that some non-relativistic particles will remain [5] even if the dynamics become non- 
linear, but in that case a full lattice field theory simulation [10, 11] is required to calculate 
the dependence of on cr*. This result could then be subsituted to Eq. (2.6) to compute 
the curvature perturbation. Simulations of this type are, however, beyond the scope of this 
paper, and therefore we will use Eq. (2.7) to obtain a rough estimate of the perturbations. 

We also assume that there is some mechanism for a perturbative decay of the curvaton. 
This is likely to be suppressed because there are no direct 3-point couplings to the Standard 
Model. There is, however, a strict lower bound on the curvaton's final effective decay rate 
(r > 10^^^ GeV) from requiring the curvaton to decay before BBN. To avoid large isocur- 
vature perturbations, the curvaton should also decay before dark matter freezes out, but this 
constraint is strongly model-dependent^ and we do not impose it. Both the resonance and 
thermal interactions are likely to occur too early to produce sufficient (. However, provided 
that some non-relativistic curvaton remains after these processes, the curvaton's energy den- 
sity can continue to grow, before perturbatively decaying and producing the majority of the 
curvature perturbation. 

Either the interactions of an initial thermal bath with the condensate, or the resonant 
production of photons by the condensate will leave a distribution of curvaton particles that 
may not be in chemical equilibrium.^ At this point, the evolution depends on the tempera- 
ture. If T <C m, then the curvaton's energy density continues to evolve as oc a~^, and a 
late perturbative decay will eventually give an adiabatic (, which could be large enough to 
match observations. If T ^ m but chemical equilibrium is never reached, then the curvaton 
energy density scales like po- oc until the temperature drops and the curvatons freeze out. 
From this point, the evolution is the same as the T <C m case, and a large C is possible. 
Alternatively, if chemical equilibrium is reached, then the number density of the curvatons 
is in equilibrium, and ( becomes adiabatic. In this case, ( is unaffected by the freeze-out or 
decay of the curvatons, and is likely to be very small (because chemical equilibrium occurs 
at an early time). 

3 Field and Charge Fluctuations 
3.1 Field fluctuations 

The curvaton is a light field, and will therefore gain perturbations during inflation, in the 
same way as the inflaton field. Assuming for simplicity that m = 0, the linearised equation 

"For reference, WIMP dark matter that decouples at T ~ 0(10) GeV gives F > 10"^'^ GeV. 
^We are careful to distinguish between kinetic equilibrium (distribution of momenta) and chemical equi- 
librium (number densities of particles). 
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of motion for a mode with comoving momentum k is 



fr2 

a(k) + 3H&(k) + ^a(k) = 0, (3.1) 

where k = |k|. The mode is initially in its vacuum state, and when it crosses the horizon, 
i.e., k/a = H, it freezes out, leading to a scale-invariant spectrum of fluctuations 

(at(k)a(q)) = (2^)3<5(k + q)Goo(A;), (3.2) 

with 

Goo(A;) = (3.3) 

For comparison with observations, we split the field into a background part which is 
the average field value in our currently observable universe, and fluctuations 6a^, which are 
local deviations from this background value. That is, within the observable universe we have 

(j(x) = (J* + (5cj=k(x). (3.4) 

The precise value of (t=k is not determined by the theory, but its typical value can be estimated 
by calculating the variance 

(Kl) = y^^^ (^Goo(^) = ^ln^, (3.5) 

where Acmb is the comoving horizon size at the time when the largest observable scales left 
the horizon, and Air is the comoving horizon size at the start of inflation. This can also be 
expressed as 

= (iVtot-iVcMB)^, (3.6) 



47r2' 

where A'tot is the total number of e-foldings of inflation, and A'^cmb ~ 60 is the number of 
e-foldings after the largest observable scales left the horizon. Although A'tot is not determined 
by the theory, we would typically expect A'tot — A'cmb ^ and therefore cr* > -ff*. When 
the non-zero mass of the curvaton is taken into account, this estimate eventually breaks down 
when the probability distribution function for o"* approaches the equilibrium solution, with 
(|(T=fp) ^ ("if) Hi [12]. For large m/H^, the approach to this equilibrium solution is 
particularly fast [13]. The initial conditions for the curvaton model are discussed further in 
[13]. 

In the same way, the variance of the field fluctuations 5a^ relative to the background 
value is 

3.2 Charge fluctuations 

One obvious question about the electrically charged curvaton is whether it generates a non- 
zero charge density with observable consequences, either during or after inflation. To address 
this, we have to calculate the charge density and charge distribution of the curvaton field. 
Because the electric charge is conserved, any charge density on large scales survives at least 
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until the relevant scale enters the horizon. Therefore, there would still be a non-zero electric 
charge density in the observable Universe due to the fluctuations of the curvaton field, and 
we have to check if it is compatible with observations. The current bound on a homogeneous 
electric charge density per baryon is qe-p < 10~^^e [14]. If the charge distribution is not 
symmetric, this bound relaxes to Qe-p < 5 x W~^^e [14]. 

A large field value does not necessarily imply high charge density. Instead, the charge 
density per physical volume of a complex scalar field a is given by 

Q = g'lm (dV) . (3.8) 

Because this involves the time derivative of the fields, it is suppressed on superhorizon scales 
where the field evolution is described the by overdamped equation 

3Ha(k) + ^aCk) = 0. (3.9) 

To calculate the statistical properties of the charge density, we write the complex cur- 
vaton field (T=K field in terms of its real and imaginary parts a = (cti + ia2)/V^, which both 
have the same correlation function 

(ai(k)ai(q)) = (a2(k)a2(q)) = (27r)3<5(k - q)Goo(fe). (3.10) 
In terms of these fields, the charge density is 

Q = ^ ((jlCJ2 - (Ti(T2) . (3.11) 

We can use Eq. (3.9) to calculate the correlation functions for the time derivatives, 
(ai(k)ai(q)) = (a2(k)a2(q)) = (2^)35(k - q)Goi(A;), 

(di(k)<>l(q)) = (a2(k)d2(q)) = (2^)35(k - q)Gii(A;), (3.12) 



and we find 



On average the charge density is zero because fiuctuations of either sign are equally likely. 
However, there will be local charge fluctuations with the two-point correlation function 

(Q(x)Q(y)) = g'^ [Gn(|x - y|)Goo(|x - y|) - Goi(|x - y|)2] . (3.14) 
Transforming Eq. (3.14) to momentum space, we obtain 

(g(k)Q(q)) =5''(2^)35(k + q) I ^[Gn(p)Goo(k-p)-Goi(p)Goi(k-p)] 

_ ^3x., A f d^P ((k-p)2-p'^' 



{27ryS{k + ci) I (3.15) 



72a4 ^ ^ ^ ' J (27r)3 p^k-p 



We can approximate this integral by noting that the dominant contribution comes from the 
two singularities p = and p = k. The integrand is symmetric, so the two poles are identical. 
Thus, expanding around p = 0, we find 

(Q(k)Q(q)) - |^(2vr)^5(k + q) y (3.16) 
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The UV divergence is cut off by k, and assuming an infrared cutoff A ^ A;, we liave 

r/'^H'^k f'^ fin n'"^ H"^ k 

(Q(kMq)) « f-^(2vr)3<^(k + q) / ^ = (2nf6{k + ci)^-^^klog- 

= {27Tf6{k + ci)GQ{k), (3.17) 

wliere we have defined 

This correlation function describes the charge density on superhorizon scales. By Fourier 
transforming it, we can write down the correlation function in coordinate space 

(Qmi^n - / = fgj in (f ) (3^19) 

where r is comoving distance. This shows that there are long-range power-law charge fluc- 
tuations whose amplitude is growing as the logarithm of the scale factor. 

In order to see whether the electromagnetic interactions have a significant effect on 
the field dynamics, we estimate the typical electric field E on the scale of the horizon. From 
Gauss's law we find that, in momentum space, the electric field due to the charge distribution 
is given by 

E(k) = (3.20) 

and therefore the variance of the electric field on the horizon scale is 

To check whether this field infiuences the dynamics significantly, we compare it to the accel- 
eration due to the vacuum energy during inflation, which is 

f = H^x. (3.22) 
Estimating ln(aH/A) ~ 100, the typical acceleration due to the electric fleld, 

Ifl = ~ 0.01 — , (3.23) 

m m 

is significant at distances \x\ < 0.01/m. In this model, we find that 10^^//* "^m < (see 
Section 6), so this effect could indeed be relevant on super-horizon scales. 

However, as long as H > nie/g' , where me is the electron mass, the electric field (3.21) is 
strong enough to create electron-positron pairs through the Schwinger process [15]. Therefore 
the charge fluctuations of the curvaton field will be screened by an opposite fluctuation in 
the electron-positron charge density. The remaining electric flelds \E\ ~ rn1/g' are far too 
weak to give rise to significant acceleration, and therefore we conclude that the backreaction 
to the curvaton evolution is likely to be negligible, meaning that the curvaton fluctuations 
have the usual nearly scale-invariant spectrum. For the remainder of the paper, this effect is 
ignored. As a result of the screening, the Universe will also be charge neutral, and therefore 
there is no conflict with the observational bounds on the charge density of the Universe. 
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4 The Effective Potential 



The large value of the curvaton gauge coupling g' gives rise to corrections to the potential, 
which can have a substantial impact on the parameter space of the model [16]. Assuming 
that the one- loop corrections are dominated by the gauge coupling 5', the effective potential 

Kff(a) = mVP + |!^H^lnf^V (4.1) 



647r2 ' ' 

where we have absorbed the self-coupling constant A into the definition of the renormalisation 
scale \JL. Therefore, the free parameters are m and /i. 

The shape of this potential is shown in Fig. 1. For small /i, the U(l) symmetric vacuum 
o" = is the only minimum. At larger ^, a second minimum with o" 7^ appears, and when 
[i is large enough, this symmetry-breaking minimum becomes the true vacuum. 

In order for a to act as a curvaton field, the effective potential has to satisfy certain 
conditions: 

1. Vacuum stability: The symmetric vacuum ex = has to be the true vacuum. Oth- 
erwise the Universe would tunnel into the vacuum with o" 7^ 0, which would break the 
U(l) symmetry spontaneously and make the photon massive. This is not allowed by 
the experimental constraints on the photon mass and gives the constraint 

9 647r2exp(l) 9 , 
^ < 3/ • (4.2) 

2. Shallow potential: In order to gain a nearly scale-invariant spectrum of perturbations 
during inflation, the curvaton field has to be light compared to the Hubble rate K^. 
This means that its effective mass meg , defined as the second derivative of the effective 
potential at the relevant field value, has to be less than the Hubble rate during inflation, 
so 

<fL = « iil- (4.3) 




Figure 1: The effective potential Kff(|cr|) of the curvaton field. The dashed line shows the quadratic 
tree-level potential, and the three solid lines show the effective potential for increasing /i, from left to right. 



-8- 



In addition, the curvaton should be subdominant compared with the inflaton potential 
because otherwise it is an inflaton. This means that the effective potential has to satisfy 



(4.4) 



Both of these conditions, (4.3) and (4.4) restrict the height of the effective potential, 
and also the maximum value of the curvaton field a*. 

3. Linearity: The curvaton has to evolve linearly both during and after inflation. During 
inflation, we have to make sure that the effective mass (4.3) is light over the whole range 
of field values present in the observable universe. This imposes the constraint 



where, according to Eq. (3.7), Sa^, ~ \/NH^. 

After inflation, we require that the potential is dominated by the mass term m^|cjp, so 



This implies that the curvaton has the equation of state of matter after inflation, making 
it easy to study the dynamics. The condition also guarantees that if a metastable second 
minimum exists, the field is on the left side of the barrier and starts to oscillate around 
the symmetric minimum. We note that it may be possible to relax the constraint (4.6), 
and have a curvaton model where the potential is initially dominated by the correction 
term. This has been discussed in Ref. [16]. 

The effect of these constraints is shown in Fig. 2 (allowed region is shaded). They are 
highly non-trivial, both on the value of , and on the mass of the field (in generic curvaton 
models, the curvaton mass m can be as low as 10~^^H^:). The combination of these constraints 
favours large curvaton mass m > 10~^H^, relatively (but not unnaturally) low field values 
o"* < O{W0)H^, and /i < 0(100)1/*. The maximum a* for various constant values of n is 
also shown, to demonstrate that larger values of cr* are only possible for larger values of /i. 

5 Evolution after inflation 

After the end of inflation, the curvaton field consists of a homogeneous condensate that 
oscillates in its potential according to Eq. (2.2). The evolution of this condensate is affected 
by its interactions with other fields, which cause it to decay into curvaton particles and other 
degrees of freedom. By curvaton particles, we mean any curvaton field modes with non-zero 
momentum, and for simplicity we assume that they have a thermal spectrum, so that their 
behaviour is completely parameterised by their number density and temperature. 

The U(l) charge of the curvaton field does not allow a direct Yukawa coupling to Stan- 
dard Model fermions, but we assume that there is some indirect decay channel, through 
non-renormalisable interactions or some beyond-the-Standard-Model fields. This allows the 
curvaton to decay perturbatively into fermions at a slow rate F which we treat as a free 
parameter. At earlier times the curvaton's interactions are dominated by its U(l) gauge cou- 
pling, which has two effects: it allows the curvaton condensate to decay non-perturbatively 



(4.5) 




(4.6) 
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Figure 2: Allowed region (shaded) for a viable curvaton model satisfying conditions 1.-3. in the text. Note 
that these conditions only guarantee the existence of a curvaton, but not that the generated perturbations 
are compatible with observations. Also shown are the maximum a* for various constant /i; vacuum instability 
rules out large n and small m. The parameter space is independent of ff, . This is because H, only determines 
the subdominance of the curvaton during inflation, which is naturally satisfied provided the other conditions 
are met. 

to photons through a parametric resonance; and thermal photons (if present) interact with 
the curvaton condensate, turning it into curvaton particles. 

These processes affect the behaviour and energy density of the curvaton field and there- 
fore the curvature perturbation through Eq. (2.7). The curvature perturbation becomes 
adiabatic when the curvaton either decays or obtains full (chemical and kinetic) thermal 
equilibrium with the dominant background. To obtain a sufficiently high amplitude of per- 
turbations, it is usually necessary for the curvature perturbation to become adiabatic very 
late. Note that the mechanism of decay can also affect the non-Gaussianity of the model. 

From Fig. 2 we can see that a^, < lOOi/*. Together with Eq. (2.7), this implies that 
in order to generate the observed amplitude of perturbations the energy fraction in the 
curvaton field r^ec must less than 10~^. Therefore we can safely assume that the curvaton's 
contribution to the energy density is subdominant and the background energy density scales 
like radiation. As discussed in Section 3.1, field values well below the Hubble rate are 
unnatural, and therefore we also set a lower limit a^, > 0.1/7*. 

We now discuss these processes in detail, beginning with the effect of a thermal bath 
created by the infiaton. 

5.1 Interactions with the thermal bath 

Because the curvaton-photon coupling g' is relatively large, the interaction of the curvaton 
with any thermal bath containing photons is significant. In particular, interactions with pho- 
tons kick curvaton particles out of the condensate, creating a thermal spectrum of curvaton 
particles. This process has the rate [18] 

Fth « Omg'^T, (5.1) 
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where T is the temperature of the thermal bath. 

This thermal interaction only converts the condensate into curvaton particles, but does 
not cause the curvaton particles to decay. However, the non-perturbative decay of the cur- 
vaton into photons (through parametric resonance) is then prevented, because this process 
requires a homogeneous condensate. If the temperature is low (T <C m), then the produced 
particles are non-relativistic and the equation of state is unchanged compared to that of the 
oscillating condensate. If the temperature is high enough (T S> m), the produced curvaton 
particles are ultrarelativistic and have the equation of state of radiation. In this case, the 
evolution depends on whether both chemical and kinetic equilibrium are reached. 

In our scenario, a thermal bath of radiation can be generated immediately after inflation 
by the decaying inflaton field, or later by annihilating curvaton particles. Let us first focus 
on the former case. Assuming that the inflaton has decayed completely to thermal radiation, 
the temperature of this radiation is given by (assuming g^, = 100) 

T ~ 0.5^/MpH. (5.2) 
According to Eq. (5.1), this gives an effective thermalisation rate 

Tth - W-'\I^H. (5.3) 

The process is efficient if the thermalisation timescale is less than the age of the Universe 
i.e., l/Fth < l/-ff. For all reasonable values of parameters, thermalisation occurs effectively 
instantaneously (with T ^ m), destroying the curvaton condensate. If the curvaton also 
reaches chemical equilibrium, this makes the produced curvature perturbation too small to 
be compatible with observations. 

Chemical equilibrium requires annihilation or pair creation of curvaton-anticurvaton 
pairs. We estimate the rate of these processes as 

rchcm ~ ]^g^2 '^th (5.4) 

which is suppressed relative to Eq. (5.3) by two powers of the coupling g' . Chemical equi- 
librium is achieved if l/Fchem *C 1/i?, which is quickly satisfied for all reasonable values of 
parameters, compared to the timescale needed to produce sufficiently large C,. For large H^, 
and m, the timescale may be comparable to that of the non-perturbative effects, discussed 
below. 

Thus, in order to have a viable model we must impose the requirement that there are 
(almost) no photons in the thermal background after inflation. This could occur either (i) 
if the inflaton decays to a hidden sector, sufficiently decoupled from the Standard Model, or 
(ii) if the inflaton is blocked from decaying until a sufficiently late time.^ In this second case, 
we require the inflaton to oscillate in a 0^ potential in order that r can grow sufficiently. We 
assume that one of these two possibilities can occur. 

5.2 Non-perturbative decay of the curvaton 

Provided we avoid thermalisation by the inflaton's decay products, the first interaction of the 
curvaton is non-perturbative production of photons, through a parametric resonance [19-21]. 

®Even a small fraction of the inflaton decaying into Standard Model particles could be enough to cause 
thermalisation of the curvaton condensate. The exact limits depend on the parameters. 
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This is particularly important because of the large photon-curvaton coupling. Note, however, 
that charge conservation prevents perturbative decay through this coupling. 

In the temporal Coulomb gauge, the equation of motion for the gauge field A in the 
homogeneous background (2.2) is 

A(i, k) + H{t)A{t, k) + + 2g''a{tf^ A{t, k) = 0, (5.5) 

where k is the comoving momentum. The gauge field has been normalised in such a way 
that for g' = it oscillates with a constant amplitude. 

Rescaling the gauge field A(t,k) with B(t,k) = a(t)^/^A(f, k) and substituting for a 
using (Eq. (2.2)) gives a Mathieu equation with time-dependent parameters, 

B"(z,k) + (Sfe(z) + 2g(z)cos2z)B(z,k) = 0, (5.6) 

where z = mt — Svr/S is a dimensionless time parameter and ' denotes differentiation with 
respect to z. The parameters are 

2i/2r(5/4)2 5'2 V 37t\-^/^ ^ g'-^al 
q{z) = ^-A^^lz+ \ -0.37^^^ 5.7 



and 



2mH^{z + 3Tr/8) 16(z + 37r/8)2 

Depending on the parameters S and q, the solutions of the Mathieu equation are either 
oscillatory or exponentially growing. The growing solutions correspond to rapid, resonant 
energy transfer from the curvaton to the U(l) gauge field. The parameter values for which 
this happens form instability bands, as shown in Fig. 3 (shaded regions are stable). In our 
case, both parameters q{z) and T,k{z) decrease with time, so that if they are initially large, 
they move through instability bands, until they leave the last instability band when g ~ 1, 
and the resonance ends. The resonance can also end if backreaction becomes so significant 
that the linear approximation (5.5) fails. The trajectory for fc = modes is shown in Fig. 3. 
Both the speed at which the solution moves to small values of the parameters, and the initial 
value depend on m and (7^,. The trajectory for modes with /c > is shallower, and thus 
modes with large k do not spend enough time in the instability bands in order to produce a 
resonance. 

In order to estimate the amount of energy transferred from the curvaton fields, we solve 
Eq. (5.6) numerically, finding the amplification factor a{z) = Ai{z, k) /Ai(0, k). This is shown 
in Fig. 4 as a function of the time variable z for two sets of parameters. If the amplification 
factor becomes large, then backreaction can no longer be ignored. We estimate that this 
happens when 
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Figure 3: Instability chart of the Mathieu equation. Shaded regions show the stable bands; white regions 
show the resonance bands with exponentially growing solutions. The solid line shows E — 2q. For fe = 0, the 
solution moves towards the origin following this line very closely. The starting position and speed at which it 
moves depend on m and cr*. Modes with k > follow a similar evolution, but along a shallower line. Thus, 
modes with higher k spend less time in the instability bands, leading to a weaker resonance. 



To approximate this, we assume that modes with k < fcres ~ \/ g'raa^f are aU amphfied by a 
factor a, and we find 



so that the resonance is non-hnear if 



^2 



(J if I n if 



(5.10) 



(5.11) 



In many cases, only a reasonably modest amplification factor a is required in order for 
the resonance to become non-linear. If the initial value q^^ w {g' aif/m)^ /2 of the resonance 
parameter q is large, the parameters q and S move slowly through a large number of resonance 
bands. Therefore the amplification factor a becomes exponentially large, and we expect that 
the dynamics become non-linear. In this case a full numerical study is necessary to determine 
the dynamics. If, on the other hand, ~ 1, then the system moves quickly through the 
resonance bands, and we do not expect significant non-linear effects. Fig. 5 shows this initial 
value of the resonance parameter. 

The value of z when the resonance either ends or becomes non-linear is z^es ^ 1000. 
It has been shown that for a non-linear resonance, non-relativistic curvaton particles are 
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Figure 4: Amplification factor a as a function of z for fc = 0, ct, = 8i?* and m = O.li/* (black, upper curve), 
m = O.JHt (grey, lower curve). The dashed line shows the nonlinearity condition (5.11). For m — O.IH, a 
huge amplification is seen and the resonance is clearly nonlinear. For m — , the amplification is much 

less dramatic and the evolution is expected to remain linear. 




10 ^ 10 ^ 10 ^ 10 

m/H. 



Figure 5: Initial resonance parameter g» in the allowed parameter space (shaded). The parameter is largest 
for large cr«, small m. A stronger, non-linear resonance is expected for larger . 



likely to remain after the resonance has completed [5]. As a very simple estimate we assume 
that a fraction / = 0.5 of the condensate's energy density is transformed to relativistic 
photons, and half remains as either condensate or non- relativistic curvaton particles. This 
assumption allows us to use the standard expression (2.7) for the curvature perturbation. 
However, simulations in a scalar model (see the bottom panels in Fig. 3 of Ref. [5]) show that 
the fraction / is actually highly dependent on the curvaton field value, and as discussed in 
Section 2, this dependence modifies the predicted curvature perturbation significantly. For 
fully reliable predictions, a non-linear field theory calculation is therefore required. 

The fraction of the original curvaton which does not decay resonantly can both interact 
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with the newly-produced thermal bath and decay perturbatively (to fermions). The final 
will include contributions from each process. There is also a direct contribution to the 
curvature perturbation generated during the resonance. Its magnitude is determined by 

r(tres) = 0.74M(^^. (5.12) 

Even by using maximum values / = 1, fj* = lOOif* and z^es = 1000, we obtain r-(tres) ~ 
10-^2(j^^/lQio GeV)2, which gives Cres ~ 10"^^, far below the observed value of 10"^. Thus, 
we conclude that the direct contribution to the observed ( from the resonant decay is negli- 
gible. 

5.3 Subsequent evolution 

Once the resonance has begun, a thermal bath of photons is created. The remaining conden- 
sate can then interact with this thermal bath. The temperature is given by 

T ^ OApll' 

^OAf/^{mt,,,y/'^ (5.13) 
and the effective rate of thermal interactions (using Eq. (5.1)) is 

/2 .1/4 

(mi 

^l/8 1^ 



rth(t > ires) = O.OlV'/^/* (mtrcs)'/' y y . (5.14) 

The thermal interactions will eventually become important when l/-ff — l/H^cs ~ l/rthj 
giving 

Evolution after the thermal interactions strongly depends on T/m and we now discuss the 
two limiting cases. 

If T <C m, then the condensate decays into non-relativistic curvaton particles. Because 
the energy density is still proportional to 1/a^, this does not affect the curvature per- 
turbation C,. Therefore, it is determined by the perturbative decay, which takes place at 
iport ~ l/r, where T is the perturbative decay rate. The curvature perturbation C is then 
given by Eq. (2.7), using r(tpert)- 

If T m, the thermalised curvatons are ultrarelativistic and the evolution depends on 
whether chemical equilibrium is reached. This is determined by tchem = 1/rchem (Eq. (5.4)), 
which should be compared to both the timescale of the expansion (l/H) and the time until 
the relativistic curvaton freezes out again (ffreeze)- Freeze out will occur when T ~ m, i.e. 

_ 0.14/V^(mt,es)^/V. 

'-freeze — 9 • lO.lDJ 

Thus, chemical equilibrium occurs if 

*chem ^ ^th and tchom ^ ^freeze " *th- (5-17) 
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In the case that chemical equihbrium occurs, the curvature perturbation ^ becomes 
adiabatic, and therefore its final value is determined by the curvaton energy fraction at the 
time of equilibration rth- Similarly to Eq. (5.12), we obtain 

Using maximum values / = 0.1, a^: = 100//*, m = H^, and Zrcs = 100, we find tth ~ 10^ and 
thus that rth ~ 10~^^(i/*/10^'^ GeV)^. This gives Cth ~ 10~^^, again far below the observed 
value of 10~^. Thus, we conclude that if the curvaton reaches chemical equilibrium with the 
thermal bath, then the contribution to C, is negligible. 

If, on the other hand, the curvaton does not reach chemical equilibrium, it will eventually 
become non-relativistic again. The ultra-relativistic period reduces the amplitude of the 
curvature perturbation C,. In order to obtain the observed amplitude, the perturbative decay 
must then be delayed compared with the non-relativistic case. In this case, is given by 

p<' = ^-^^-fil^u — • ^^^^^ 

t ' f-rreeze 

The calculation of Q then follows exactly the same procedure as for the non-relativistic case. 
The only difference is that a somewhat smaller F will be required to obtain sufficient and 
in some cases this value of T could be ruled out by the BBN limit. 

6 Observational constraints 

We explore the parameter space numerically, focusing on the observables C (curvature per- 
turbation), /nl (non-Gaussianity) and n (CMB spectral index) in turn. 

6.1 Curvature perturbation 

In general, the predicted amplitude (2.6) of the curvature perturbation depends on the 
perturbative decay rate F, because it sets the value of tdec- Where possible, we fix T in order 
to obtain C, = 10~^, however there is a constraint on T. The curvaton must decay before 
BBN, which means that extremely small values of T are not allowed. This restricts the low 
-ff*, low m region of parameter space, because r oc (cr*/ff*)^(i7*/Mp)^(m/r)^/^. Thus low 
m, low i/* gives low r. For < 10^ GeV the requirement C, = 10~^ means that there is 
no available parameter space left. For 10^ GeV < 10^ GeV, the parameter space is 

reduced (Fig. 6). For H^, > 10^ GeV, there is no effect on the parameter space. 

6.2 Non-Gaussianity 

The observational limits on the non-Gaussianity of the curvature perturbation impose further 
constraints on the parameter space of our model. To obtain a rough estimate, we ignore the 
effect of the parametric resonance and use the standard expression (2.8) for /nl- In the 
parameter space with the correct C, we can use Eq. (2.7) to find r^eo giving 

rdec = 37rC-^. (6.1) 

Using cj*/i/* < 100 gives, 

/nl > 130. (6.2) 
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Figure 6: Allowed region for a viable curvaton model which produces — 10^^. For Ht. > 3 x 10^, the 
requirement is satisfied in the entire parameter space (all three shades of grey). The size of the allowed region 
reduces as Ht reduces. For Ht = 10^ GeV, the allowed region comprises of the two darkest grey regions; for 
i/t = 2 X 10* GeV it is only the dark grey region. For Ht < 10* GeV there is no allowed parameter space. 



The 95% confidence level WMAP exclusion is -10 < /nl < 74 [22]. Being within this 
range would require o"^, > 180H^, so this would appear to rule out the model. However, 
this conclusion is not warranted because, as discussed in Section 2, the parametric resonance 
modifies the predictions significantly. Thus, numerical lattice simulations are necessary. 

6.3 Spectral index 

Another observable parameter that could be used to constrain the model is the spectral index 
n. Unlike C and /nLj it depends on the specific model of inflation, and therefore the results 
are less generally applicable. As an example, we consider the simple monomial potential 

V{cl)) = X^cl)\ (6.3) 

This is an obvious choice in our case, because we require a non-thermal inflaton back- 
ground after inflation that behaves like radiation. Requiring A'^ e-foldings of inflation gives 



/>2 ~ 8NMp. The spectral index for the curvature perturbation is 



n = 1 - 2einf + 2r/curv 

2 2m2 , ^ 

The current observational limits (two-sigma) are 0.939 < n < 0.995 [22]. For A^ = 50 and 
A^ = 60, n is within the WMAP limits provided m < 0.2i?^,. This would rule out the large m 
region of the parameter space, but we stress that this is dependent on the model of inflation. 

6.4 Cosmic strings 

An additional interesting physical effect could occur because the curvaton field effectively 
has a non-zero value within any Hubble volume, therefore breaking the Standard Model 
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U(l) symmetry spontaneously. For topological reasons there will be curves in space where 
the curvaton field vanishes, in very much the same way as in a cosmic string. At the end 
of inflation, the curvaton fields starts to oscillate, and these strings dissolve, but because 
the field value is zero at the string locations, they can potentially influence the curvature 
perturbation on cosmological scales. Furthermore, just like ordinary cosmic strings, these 
strings carry one quantum of magnetic flux $o = 27r/(7', which can also have an effect on 
the reheating dynamics locally. An investigation of these effects is currently underway, using 
lattice simulations. 

7 Conclusions 

We have demonstrated that it is possible to have a consistent model of the early Universe in 
which a scalar field charged under the Standard Model U{1) weak hypercharge plays the role 
of the curvaton, generating a nearly scale-invariant spectrum of curvature perturbations with 
the observed amplitude. Besides curvature perturbations, the curvaton charge gives several 
other potentially observable effects. In principle, field fluctuations during inflation generate 
significant electric charge fluctuations on superhorizon scales, but we found that these are 
suppressed by the screening provided by Schwinger pair creation of other charged particles 
such as electron-positron pairs during inflation. 

The standard calculation for this model predicts relatively a high Hubble rate during 
inflation, > 10® GeV, and significant non-Gaussianity (/nl ^ 130), which would rule it 
out at 95% confidence level by WMAP. However, this ignores a period of parametric resonance 
between the curvaton and photon fields, whose effect on the curvature perturbation can only 
be calculated with numerical lattice field simulations. There is also freedom to relax some of 
the constraints in Section 4, such as allowing a meta-stable vacuum or permitting a quartic 
term to initially dominate the curvaton potential. 

Infiaton dynamics after inflation are also strongly constrained to avoid thermal photons 
destroying the curvaton condensate. In practice, this means that inflaton reheating should 
be substantially delayed, and that the inflaton should oscillate in a quartic potential after 
inflation (the inflaton should behave like radiation so that the relative magnitude of the 
curvaton can grow). 

The existence of an extra U(l) charged fleld, with a non-zero value, could also play a 
significant role in the electroweak phase transition. However, investigating these possibilities 
is beyond the scope of this paper, and we simply raise them here as potential directions for 
future research. 

In summary, this paper presents an interesting model that concretely links the dynamics 
of the Standard Model of particle physics with the very early Universe. In the era of new 
data from both LHC and Planck, the mechanism presented in this paper deserves further 
investigation. 
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